Tracking quasi-classical chaos in ultracold boson gases 
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We study the dynamics of a ultra-cold boson gas in a lattice submitted to a constant force. We 
track the route of the system towards chaos created by the many-body-induced nonlinearity and 
show that relevant information can be extracted from an experimentally accessible quantity, the 
gas mean position. The threshold nonlinearity for the appearance of chaotic behavior is deduced 
from KAM arguments and agrees with the value obtained by calculating the associated Lyapunov 
exponent. 
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Recent advances in the physics of cold atoms paved the 
way for the investigation of fundamental quantum prob- 
lems with unprecedented cleanness. "Quantum chaos" is 
one of the most fascinating among these problems, be- 
cause in such case the correspondence between the dy- 
namics of the quantum system and its classical coun- 
terpart is nontrivial. The reasons for that are twofold: 
First, quantum particles obey Heisenberg inequalities, 
hence their dynamics cannot be described in terms of 
phase-space trajectories. Second, sensitivity to initial 
conditions is not observed, as the Schrodinger equation 
is linear. For these reasons, quantum chaos has often 
been defined as "the behavior of a quantum system whose 
classical limit is chaotic". Whereas this is a reasonable 
definition, it is clear that the actual quantum dynamics 
has no direct relation to classical chaos. Much work thus 
concentrates in finding "reminiscences" of classical chaos 
that might survive in the quantum system, the so-called 
"signatures" of quantum chaos [H . 

The realization in 1995 of the first Bose-Einstein con- 
densate with laser-cooled atoms 0, 0, 01 opened a new 
way for investigating (truly) nonlinear dynamics in quan- 
tum systems. In an "ideal" (i.e. in the zero temperature 
limit) Bose-Einstein condensate (BEC), the atoms are 
indistinguishable and form a mesoscopic object which 
can be described by a "collective" single wavefunction. 
The EEC's dynamical behavior is then described by 
the Gross-Pitaevskii equation which includes a nonlin- 
ear term due to atomic interactions The solutions of 
such equation can - and do - present sensitivity to initial 
conditions, leading to "classical-like" instabilities, a pos- 
sibilitythat attracted much attention from both theoreti- 
cians M 0, S, @, 11 , H, 13, 3 and experimentalists 

Mm- 

In order to make these ideas clear, let us first consider 
a simpler system, the kicked rotor, whose experimental 
realization with laser cooled atoms [13 lead to an im- 
pressive burst of experimental work [isl . 19l . 2ol . 21 , [i^, 
Q, [H, [IHl • This system is formed by a particle periodi- 
cally "kicked" by a sinusoidal force. The classical dynam- 
ics is characterized by a single parameter K, which is 
the normaHzed amplitude of the potential whereas 



the "quanticity" of the quantum system is described by 
a "normalized Planck constant" k (the dynamics become 
classical as ^ — > 0). 

In order to characterize the dynamics, one may chose, 
for instance, the particle's square momentum as the rel- 
evant quantity (p^(i) in the classical case, (p'^it)) in the 
quantum case) . By taking the Fourier transform of the 
above quantity and pinpointing the frequencies appear- 
ing in the evolution of the system for different values of 
K and fixed initial conditions, we can characterize in a 
common language the classical and the quantum versions. 
The result is plotted in Fig. [TJ 

The classical dynamics. Fig. [T^, shows resonances (or 
frequency-lock events) each time a frequency become very 
close to another one. For K > 0.8, the frequencies in the 
spectrum "coalesce" to form of a dense spectrum, which is 
a signature of a chaotic dynamics. In the quantum case, 
Fig. [TJd, the spectrum is richer and present numerous 
frequency crossings, but no dense spectrum is observed. 

Let us now consider a one-dimensional BEC of parti- 
cles of mass M placed in a sinusoidal lattice formed by 
counterpropagating laser beams of wavelength Al, and 
subjected to a linear force. The system is described by 
the Gross-Pitaevskii equation (GPE) 



'' dt 



2m 



+ Vo cos{2'Kx) +Fx + g j Tp{x, t), 

(1) 

with lengths measured in units of the lattice step d = 
Al/2, energy measured in units of the recoil energy Eu = 
h^kl/{2M) (with kL = 2t:/\l), time in units of h/ER, 
the force F in units of Er/d, P = —id/dx, the reduced 
mass is TO = 7r^/2 and g is the (ID) nonlinear parameter. 
The eigenstates of the linear part of the Hamiltonian ([TJ , 
obtained by settingg = 0, are the so-called Wannier- 
Stark (WS) states [10, W^- In order to simplify the 
discussion, let us we suppose that the BEC energy is low 
enough that its wavefunction -0 can be expanded only 
on the lowest-energy WS state for each potential well, 
hereafter noted as (pn(a;) (corresponding to the n*^ well). 
These states are centered at each potential well, obey 
the symmetry relation ipn+m{x)= ipn{x — m), and form a 
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Figure 1: Classical and quantum dynamics of the kicked ro- 
tor. For each value of we numerically calculate 
(classical case) and {p^{t)) (quantum case) up to 2000 kicks, 
perform the Fourier transform and pinpoint the frequencies 
whose Fourier amplitude is above a threshold (here 1/500 of 
the maximum amplitude), (a) Frequencies in the classical 
dynamics: A dense spectrum appears for K > 0.82 (initial 
conditions a; = and p = 0.56). (b) Frequencies in the quan- 
tum dynamics: The spectrum is complex, but always discrete 
{k = 2.89, the initial state is a Gaussian in momentum space 
of FWHM 5K centered at p = 0.56). 



ladder of eigenenergies £"„ = nujB, where ujb = Fd/h {F 
in our normalized units) is the Bloch frequency. Thus, 
putting (p{x) = ipo{x) 



^(a;,i) =^/^e-*«"V'(x- 



(2) 



where the eigenstates population /„ and phase 6'„ are real 
functions of the time. 

Following the approach of we insert Eq. ([2]) in 
Eq. Ilj, and obtain Hamilton equations for the popu- 
lations and phases. The Hamiltonian then appears as a 
sum of an integrable part Hq and a non-integrable per- 
turbation Hi: 



H = Ho{I)+eHi{I, 



(3) 



with e oc gxoi, where xoi = J dxtp^ {x)ip{x — i) is a mea- 
sure of the superposition of eigenfunctions correspond- 
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Figure 2: Frequencies present in the spectrum of {x{t)) vs g 
(the threshold is 1/100 of the maximum amplitude) . For small 
g, the Bohr frequencies [cf. Eq. ([4|)] are lou) = i^s — 0.4t;xoo, 
ci^o-i — LOB + 0.55(/xoo and c^i-i — 2ujb + 0.15t;xoo, corre- 
sponding to branches a, b, and c, respectively. Harmonics 
of uJio are also seen (branches d and e) and r2i,_i,o (branch 
f). For 0.22 < (? < 0.26 one observes dense-spectrum win- 
dows corresponding to quasi-classical chaos. Parameters are 
ujb ~ 0.25, V = 5 (xoo = 2); initial conditions are lo ~ 0.65, 
h = 0.25, 00 = 9-1 ^0,61= TT. 



ing to different wells. For small e, the system is quasi- 
integrable and fits into the general frame of the KAM 
theorem Therefore, the quantum-coherent evolution 
of a BEC displays classical-like KAM-structured chaos. 
We call this, "quasi-classical" chaos. 

If we temporarily neglect the nonlinear term in Eq. ([ij , 
we see that the phases evolve as 6n{t) = 9no + nujBt. 
To the first order in e, the effect of the nonlinearity is 
to introduce a population-dependent correction to this 
phase, producing a harmonic evolution with frequencies 
ujn = nujB +9XoaIn- The intensities are constants of mo- 
tion, In{t) = In{t = 0) and the "self-coupHng" parameter 
Xoo depends only Vq and F. The BEC dynamics is thus 
governed by the Bohr frequencies 



UJ„ 



uJn - uj„i = (n - m)LL>B + 9Xoo{Iri 



(4) 



In order to simplify further our description, let us restrict 
the dynamics to three adjacent potential wells, i.e. we 
set /„ = if rt ^ -1,0,1 in Eq. The system is 

then four-dimensional, the dynamical variables being two 
populations (since /_i + /o + /i = 1) and two relative 
phases. 

The onset of chaos in this system is shown in Fig. [21 
which is the equivalent of Fig. [T] for the BEC dynam- 
ics, with the difference that we used the average position 
{x{t)) instead of (p^it)"^. Strikingly, the plot resembles 
more closely to the classical (Fig. [T^) than to the quan- 
tum (Fig. [TJd) dynamics of the kicked rotor. 

We can interpret the main features of Fig. [2] from sim- 
ple arguments. At 5 = there are only three Bohr fre- 
quencies in the model: Two of them are degenerated, 
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wo-i = wio = i^B, and correspond to energy differ- 
ence between neighbor wells, whereas uj-n = 2uib corre- 
sponds to next-to-neighbor energy difference. The low-g 
structure can be understood using Eq. The frequen- 
cies present around g = 0.05, for instance, in Fig. [2] are 
of the form 

VLpqr = puJa-i + qt^io + ruji^i (5) 

with p, q, r integer. Namely, from top to bottom, we 
found wi-i, 2ijJio, cjo-i, i^w and ujq-i — loiq. These 
frequencies are only weakly perturbed by the nonlinear- 
ity, which manifests itself by the slight curvature of the 
branches. A KAM-type expansion in powers of e shows 
that the weight of a frequency flpqr in the spectrum is 
proportional to g/^pqr [28|, thus, the higher the value 
of g, the larger the number of frequencies that will be 
present (remember that Fig. [2] displays the frequencies 
whose amplitude is above a threshold), and smaller fre- 
quencies will be favored. The intersection of two branches 
correspond to a resonance, that is, the condition fJpq,. = 
is fulfilled for some p, q, r. Close to a resonance the KAM 
perturbation theory breaks down, and non-integrable be- 
haviors appear that may lead to chaos. In Fig. [2] the 
arrows indicate the resonances fliqo = (g = —2) for 
g w 0.093, and q — —3, —4, —5 for, respectively, g « 0.14, 
g = 0.17, g = 0.196 . For higher g values, resonances with 
higher and higher q values become significative, as seen 
on the left of the chaotic region: Each frequency crossing 
produces a multi-frequency structure (indicated by the 
arrows). Finally, for g > 0.2 one observes at least four 
successive windows of dense spectrum, corresponding to 
a chaotic behavior. 

Fig. [2] represents a "local route" (in the sense that its 
detailed geometry depends on the initial conditions) to 
quasi-classical chaos. The above argument suggests how- 
ever that this route is characteristic of KAM systems, and 
globally independent of initial conditions. Its universal 
character is confirmed by comparing it to Fig. [l^, where 
an analogous structure is observed in a completely dif- 
ferent KAM system, the classical kicked rotor. We have 
thus put into evidence a route to KAM chaos in a non- 
linear quantum system, which is potentially observable 
with state-of-art experiments. 

In order to confirm our conclusions, we have also cal- 
culated the maximum Lyapunov exponent (MLE) as- 
sociated to the dynamics, which is a direct signature 
of the sensitivity to the initial conditions, and thus of 
chaos. In order to calculate MLEs in our quantum sys- 
tem, we adapted the classical Jacobian method [i^. We 
represent the system evolution by a trajectory in a six- 
dimensional "generalized quantum phase-space", formed 
by the real and imaginary parts of each WS state co- 
efficient, i.e. ^/T^cosOn and v^sin0„, for n = —1,0,1 
(this is numerically more stable than using /„ and 0„), 
and calculate the divergence of neighbor WS states, from 
which we can extract the MLE. The result is presented in 




Figure 3: Maximum Lyapunov exponent (in units of Tg^ — 
uib/2-k) as a function of g. One observes a good correspon- 
dence between the existence of a non-zero Lyapunov exponent 
and the presence of dense-spectrum regions in Fig. [2] The 
observed dense-spectrum is thus a reliable signature of the 
chaotic behavior. Same parameters and initial conditions as 
in Fig. El 



Fig.m One sees a clear transition to chaos for g ~ 0.194, 
which is in good agreement with the value that can be 
deduced from Fig. [2l We can also identify the four suc- 
cessive chaotic windows. 

For given initial conditions, we can obtain an estimate 
for the critical value of g at which chaos appears. KAM 
theory shows that chaos generally appears along a sep- 
aratrix. In our system there are two main kinds of tra- 
jectories 0: "Passing" trajectories correspond to Bloch 
oscillations slightly perturbed by the nonlinearity (for the 
low values of g we are considering) ; they appear when the 
three populations are comparable. Bound, periodic or- 
bits correspond to a motion essentially confined at a po- 
tential well, and appear when one population dominate 
the others. Between these two kinds of trajectories, there 
is a separatrix. For a given trajectory (i.e. for fixed ini- 
tial conditions), chaos appears when the changing in the 
value of g brings it close to the separatrix. The condition 
for that can be simply estimated by confining the BEC 
to only two wells (/_i =0), in which case the system is 
integrable. This "reduced" system does not exhibit chaos, 
but it has a phase-space structure analogous to that of 
the 3-wells, displaying bound and passing trajectories, 
and, between them, a separatrix. In the reduced system 
we can calculate both the energy Eq corresponding to 
the trajectory (with /_i = 0) and the energy Eg of the 
unstable point to which the separatrix is connected. We 
can then numerically determine the value of g for which 
Eq = Es, which gives the critical value. Fig. [4] shows 
the result for various initial conditions (soHd fine). The 
dotted fines indicate the limits of the chaotic region as 
inferred from the Lyapunov exponent calculation in the 
3-well system. The agreement is very good, even when 
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Figure 4: Determination of the critical value of g. The solid 
lines with triangles are the estimated critical value of g, given 
by the condition that the trajectory energy equals the sep- 
aratrix energy (cf. text). The dotted lines with crosses are 
the limits of the chaotic region, where the Lyapunov exponent 
is non-zero. The parameters are the same as in Fig [21 with 
6() — 9i — ^ . 



the chaotic zone is narrow. 

In conclusion, we have characterized the phenomenon 
of quasi-classical chaos using an experimentally accessible 
signal, the mean position of the boson gas. We have stud- 
ied a "local route" to quasi-classical chaos and shown that 
the information obtained from this approach agrees well 
with that provided by a hallmark signature of chaos, the 
sensitivity to initial conditions, quantified by a positive 
Lyapunov exponent. The understanding of the structure 
of such a route allowed us to determine the critical value 
of the nonlinearity parameter, which is in good agreement 
with the one deduced from the calculation of the Lya- 
punov exponent. We think that the present work might 
stimulate an experimental observation of quasi-classical 
chaos, which would, in turn, stimulate new investigations 
on the nature of quantum chaos. 

The authors are happy to thank M. Lefranc for fruit- 
ful discussions and for his help with the calculation of 
Lyapunov exponents. 
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